A parity conserving dimer model with infinitely many absorbing states 
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We propose and study a model where, for the first time, two aspects are present: parity conser- 
vation and infinitely many absorbing states. Whereas steady-state simulations show that the static 
critical behaviour is not affected by the presence of multiple absorbing configurations, the influence 
of the initial state associated with the presence of slowly decaying memory effects is clearly displayed 
in time dependent simulations. We report results of a detailed investigation of the dependence of 
critical spreading exponents on the initial particle density. 
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Various nonequilibrium models exhibiting phase tran- 
sitions from an active to an absorbing phase have been 
extensively investigated over recent years 0-^] . 

Most of them [p|~^0[ have been shown to belon g t o 
the universality class of directed percolation (DP) 
These studies gave support to the conjecture made by 
Janssen and Grassberger 12 1^] that continuous transi- 
tions to a unique absorbing state generically fall in the 
DP class. The same applies, as far as the static critical 
behaviour is concerned, to various models with multiple 
absorbing states [ld[-|l^]. However, the dynamic critical 
properties of these systems are found to be non-universal: 
in fact, spreading critical exponents depend on the initial 
density and obey a generalized hyperscaling relation de- 
rived by Mendes et al. |l6|]. Some systems, first claimed 
to have a different critical behaviour, were later included 
in the DP class, upon a more careful analysis 

For sometime the only accepted exceptions to the DP 
universality class were the models A and B introduced 
by Grassberger and collaborators 17,lq]: they both have 



a doubly degenerated absorbing state and modulo 2 par- 
ity conservation of kinks (-00- and -11-s). Since then, a 
growing number of models whose critical behaviour falls 
in the parity conserving (PC) universality class have been 
studied. One of the first to be considered was the branch- 
ing annihilating random walk with an even number m 
of offsprings (BAW) @Jl|^|: in this case, the parity 
of the particles is conserved and there is a single ab- 
sorbing state. Another system extensively studied is the 
nonequilibrium fsing model with combined spin-flip and 
spin-exchange dynamics (NEK1M) here, the absorb- 
ing state is symmetrically doubly degenerated and the 
kinks have local parity conservation, similarly to what 
happens with the A and B models referred above. Other 
models have since been discovered and shown to display 
PC critical behaviour: the interacting monomer-dimer 
model ]23| and generalizations of the Domany-Kinzel au- 
tomaton with n equivalent absorbing states (2^] . In these 
models, parity is not strictly conserved but there is com- 
plete symmetry among absorbing states; actually, if this 
symmetry is broken, DP critical behaviour is recovered. 

Therefore, the question about the essential feature be- 
hind PC critical behaviour is not yet settled. Whereas 



in BAW models mass conservation of modulo 2 is deter- 
minant, in the other cases symmetry among absorbing 
states plays the relevant role. In the NEKIM model, par- 
ity conservation in the number of kinks is not sufficient 
to ensure PC-like behaviour, if such symmetry is broken 
II- 

Recently, Inui and Tretyakov [ p5| presented a new ver- 
sion of a contact process with parity conservation, which 
displays a phase transition for m = 2, contrarely to what 
happens with the original version of BAW . They claim, 
on the basis of Pade aproximants and also numerical sim- 
ulations, that the order parameter exponent (3 should be 

1, while the values estimated by different authors are 
consistently smaller (l9|]2|I|]. 

In this Letter, we consider a model that resembles the 
pair-contact process (PCP) introduced by Jensen Q, 
but where the number of dimers is conserved modulo 

2. Dimers cannot be generated spontaneously, therefore 
any configuration with isolated particles and empty sites 
is absorbing, and there is an infinite number of such con- 
figurations. This enables one to study the way initial 
configurations affect the spreading critical exponents in 
a system that belongs to the PC-universality class. 

Initial conditions were shown |2(| to affect time- 
dependent critical exponents that describe the damage 
spreading transition (DS) in cases where this one coin- 
cides with the ordinary critical point; in DS, one follows 
the evolution of two replicas, so this can be seen as a mul- 
ticomponent system, and any state where both replicas 
are identical is considered as absorbing. 

We start by briefly reviewing the dynamic rules for the 
PCP. This is a lattice model which can be described by 
a 2-state variabe Oi = 0,1. Nearest-neighbour pairs of 
particles (dimers) annihilate each other with probability 
p or create, with probability 1 — p, a particle at one of the 
adjacent (vacant) sites to the dimer. The annihilation 
of a pair may imply the loss of one or two other pairs, 
if one (or two) of the nearest-neighbour sites to the cho- 
sen dimer happen to be occupied; analogously when a 
particle is created adjacent to a dimer, the number of 
pairs increases by one or two, depending on the occu- 
pancy of the other nearest neighbour to the site that is 
being occupied. It is therefore clear that there is no par- 
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ity conservation in the number of dimers and the model 
displays a DP-like phase transition from an active state 
with a nonzero concentration of dimers (for p < p c ) to a 
phase with infinitely many absorbing states (for p > p c ). 
Critical spreading exponents are non-universal depend- 
ing upon the nature of the initial configuration. 

In the one-dimensional model we now introduce, a 
phase with infinitely many absorbing states is also 
present. On the other hand, dynamic rules are such that 
there is parity conservation in the number of dimers. This 
seems to be the most relevant feature; indeed, we find 
that the static critical behaviour of this model is PC- 
like, thus it is not affected by the presence of multiple 
absorbing states. 

The basic processes are again annihilation and cre- 
ation, which are attempted with probabilities p and I —p, 
respectively. If one chooses to annihilate, then one dimer 
(represented by • = •) is selected at random and one 
looks for an adjacent dimer. If the pair of dimers is 
surrounded by empty sites (o), then the two dimers 
are annihilated and the respective sites become empty 

(— o — • = • = • — o > — o — o — o — o — o — ); otherwise, 

annihilation is produced by simply vacating the site that 
is shared by the two dimers, leaving the other sites un- 
changed (— o — • = • = • = •— — > — o — • — o — • = 
•— ). In this way, the number of dimers is conserved 
modulo 2. In case of a creation attempt, then if the 
nearest and next-nearest-neighbour sites of a selected 
dimer are, respectively, vacant and occupied, a pair of 
dimers is produced by simply occupying that vacant site 
( — • = • — o — • — o— — > — • = • = • = • — o — ); 
another possibility requires the presence of three vacant 
sites adjacent to a dimer, in which case a pair of dimers 
is created by filling the nearest and next-nearest neigh- 
bour sites and leaving the third vacant site unchanged 
(—• = • — — — 0— — >—• = • = • = • — o—). In 
order that the presence of isolated particles may result 
in the enhancement of possible growth, we introduce a 
parameter a, such that the rate of creation by the latter 
process is a fraction a of the total creation rate. In the 
present work, we have fixed a = 1/2. 

With these rules, diffusion is only indirectly present. 
One can of course enlarge the parameter space to include 
diffusion, but we have not done it in the present work. 
In the absence of such processes, creation of dimers is 
inhibited in configurations where a pair of vacant sites is 
surrounded by occupied sites. Initial configurations have 
therefore been chosen such that this situation does not 
occur in the dynamical process. Also, the annihilation 
of isolated dimers is here possible only through suces- 
sive processes of creation and annihilation. This affects 
particularly (and slows down) the dynamical processes 
leading to absorbing configurations. 

The order parameter of the system is the concentration 
of dimers, which vanishes algebraically as p approaches 
the critical probability p c : 

p~(Pc-pY (i) 



where /3 is the order parameter exponent. 

Simulations were started with a fully occupied lattice 
and the concentration averaged over a long period of time 
once the steady state has been reached. 

In fig. 1 we show a log-log plot of the steady state con- 
centration as a function of p c —p for system size L = 2000; 
time varied from t = 5000 to t = 2 x 10 5 MCS closest 
to Pc(= 0.2990(10)) and we averaged over around 1000 
independent samples which had not yet entered the ab- 
sorbing state. From the slope of the data we estimate 
j3 = 0.98(5), a value which agrees with the results of 
p5| and is slightly above the values obtained by other 
authors for models in the PC universality class. Deter- 
mining the critical point by steady state simulations re- 
quires an accurate measurement of the order parameter 
near criticality, which becomes rather difficult due to the 
critical slowing down; a scaling region of over one decade 
in A = p c — p can clearly be observed in our data and 
leads to the above estimates for p c and f3. 

We have complemented this study with a finite-size 
scaling analysis based on the ansatz that the order pa- 
rameter depends on system size L through the ratio of L 
and the correlation length £ ~ A~ l/± : 

p(p,L)~L-^/(AL 1 ^x) (2 ) 

(with f(x) oc x 13 for x — > 00, so that (1) is recovered 
when L — > 00). 
Accordingly, 

p(p e ,L)ocL-^ (3) 

In Fig. 2, we show a log-log plot of p(p c , L) as a func- 
tion of L; a straight line is obtained, from the slope of 
which we estimate (3/v± — 0.54(4), to be compared with 
0.48 and 0.50 as obtained in |23| and respectively. 

We also report results of critical spreading, i.e. the 
evolution of the critical system from a configuration that 
is very close to an absorbing state. In the present model, 
there are many absorbing configurations and some criti- 
cal spreading exponents are indeed dependent upon the 
density of particles in the initial state. The quantities 
that are usually considered in these studies are the sur- 
viving probability P(t), the number of dimers n(t) av- 
eraged over all runs, and the mean-square distance of 
spreading R 2 {t) averaged over the surviving runs. At 
criticality, these quantities obey, in the long time limit, 
P(t) ~ t- s , n(t) ~ t"i and R 2 (t) ~ t z , and the corre- 
sponding exponents can be obtained from the straight 
lines shown in double-logarithmic plots of the quanti- 
ties against time. More precise estimates are usually ob- 
tained by looking at local slopes, for example —S(t) = 
ln[P(t)/P(t/m)]/ln(m). In a plot of the local slopes vs 
1/t, the critical exponents are given by the intercept of 
the curves for p c with the vertical axis, whereas curves 
for p > p c (p < p c ) veer downward (upward). This often 
constitutes a rather accurate method for the determina- 
tion of p c - It is certainly the case when one or just a few 



2 



absorbing states are present [ p3|j27| ] . Also in some sys- 
tems with infinitely many absorbing states, the critical 
point was shown not to depend on the initial configura- 
tion and can therefore be located by the procedure we 
just described In other cases, indications of a slight 
dependence of the critical point (as determined by the 
above time-dependent analysis) on the initial configura- 
tion have been found [2^p9[| and can be atributed to 
slowly decaying memory effects pf]| . This seems to occur 
also in the present case. 

In these simulations we started the system in a way 
that a sublattice of alternating sites is vacant except for 
one central site, and the sites of the other sublattice are 
occupied with probability q; the two nearest neighbours 
of the central site are also occupied, thus constituting a 
perturbation of exactly two dimers. The size of the lat- 
tice was chosen such that the spreading never hits the 
boundaries. For different q values, a number of indepen- 
dent runs, typically 10 7 , were performed, up to 4000 time 
steps each, and various values of p. 

A local slope analysis for q = 0.43 leads to the estimate 
p c = 0.2995(5), which is fully consistent with the p c esti- 
mate from steady-state simulations, above; the estimates 
for the dynamic scaling exponents are i}' = 0.005(5), 
5' = 0.291(6), z'/2 = 0.570(6) and agree well with 

§vious results for models in the PC universality class 
(ri = 0.000(1), 6 = 0.285(2), z/2 = 0.571(1)). This 
is similar to what happens in the pair-contact process: 
DP time-dependent exponents are found if one studies 
spreading of a perturbation to the "natural" configura- 
tion, the one spontaneously generated by the critical dy- 
namics. The estimate q nat = 0.43(2) (q not meaning the 
value of q that corresponds to the density in the natural 
configuration) is in good agreement, within numerical ac- 
curacy, with the one obtained by generating the natural 
configuration in samples of size L = 2000, starting with 
a full lattice and using periodic boundary conditions. 

For different q values, power laws still hold, but with 
different critical spreading exponents, and a slight shift 
of the critical point. In Fig. 3 we show the local slope 
analysis for q = 0.2, leading to the following estimates: 
6' = 0.366(6), 7/ = -0.069(5), z'/2 = 0.567(6),^ = 
0.3000(5). The generalized hyperscaling relation rj' + 8' — 
z'/2 = —6 [16| is well satisfied by these values. 

Our results for the g-dependence of p' c ,5',T)' and z'/2 
are given in Table 1. A check of the generalized hyper- 
scaling is also included. Indeed, the exponent z' does 
not present a significative dependence on the parameter 
<7, given the numerical errors. Then, as expected, the 
exponent governing the population growth in surviving 
critical trials, 8' + 77', does not depend on the initial par- 
ticle concentration. 

In conclusion, we have performed a similar study to 
what has been done before jL4| for the PCP (and other 
models with multiple absorbing states |l6| ) but now 
a system with parity conservation (in the number of 
dimers) is for the first time investigated. A field the- 
ory formalism to appropriately describe this situation is 



now required. 

Whereas the presence of multiple absorbing states does 
not affect the static behaviour, spreading critical be- 
haviour is expressed by power laws whose exponents de- 
pend on the initial particle density (= q/2). There is also 
a monotonic shift in p' cl which, when found in other sys- 
tems, has been atributed [ p9| to slowly decaying memory 
effects displayed by the non-order field. The generalized 
hyperscaling relation is verified. 

A possible mapping of this system to a SOC model |3l| ] 
is currently under study. 
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FIG. 1. Logio-logio plot of the concentration of dimers, p 
versus p c — p, with p c = 0.2990. 
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FIG. 2. Logio-logio plot of p(p c , L) versus L (L = 250, 500, 
1000 and 2000). 
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FIG. 3. Local slopes —S'(t) (upper panel), rj'(t) (mi- 
dle-panel) and z' /2(t) (lower panel) for q = 0.2. Each panel 
contains five curves with, from bottom to top, p — 0.2980, 
0.2990, 0.3000, 0.3010 and 0.3020. 
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g K 



5' g z'/2 j + 5'-z'/2 + 5 



DP - 0.1597(3) 0.312(2) 0.632(1) 0.000(1) 

PC - 0.285(2) 0.000(1) 0.571(1) -0.001(1) 

0.1 0.3005(5) 0.400(6) -0.106(5) 0.570(6) 0.01(2) 

0.2 0.3000(5) 0.366(6) -0.069(5) 0.567(6) 0.02(2) 

0.43 0.2995(5) 0.291(6) 0.005(5) 0.570(6) 0.01(2) 

0.45 0.29935(5) 0.284(6) 0.014(5) 0.570(6) 0.01(2) 

0.5 0.2990(5) 0.271(6) 0.033(5) 0.572(6) 0.02(2) 

0.6 0.2990(5) 0.239(6) 0.062(5) 0.573(6) 0.01(2) 

0.8 0.2985(5) 0.181(6) 0.126(5) 0.576(6) 0.02(2) 

0.9 0.2985(5) 0.156(6) 0.154(5) 0.579(6) 0.02(2) 



TABLE CAPTION: The table shows the g-dependence of 
the critical parameter p' c and critical exponents 6', rf and 
z'/2. A test of generalized hyperscaling is also included. 
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